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We show that different non-conventional superconductors have one fundamental 
feature in common: pair eigenstates of the Hamiltonian are repulsion-free, the W = 
pairs. In extended Hubbard models, pairing can occur for resonable parameter values. 
For (N, N) nanotubes the binding energy of the pair depends strongly on the filling 
and decreases towards a reduced but nonzero value for the graphite sheet N — > oo. 



There is experimental evidence that the critical Temperature T c in alkali-graphite intercalation compounds 
(GIC) CiM (where M is a given alkali metal) grows as x decreases [1]. Under high-pressure, high metal 
concentration samples such as CeK, C3K, C4Na, CaNa, C2Na, C2Li have been synthesized; for C2Na the value 
of T c is 5 K while for C2Li, T c =1.9 K; quite recently Potassium [2] and Lithium [3] have been intercalated 
also in single- and multi-wall carbon nanotubes [4] up to high concentration (the highest metal concentration 
was obtained with Lithium in C2Li) and a net charge transfer was observed between the alkali-metals and 
the carbon atoms. The alkali-metals cause little structural deformation, but increase the filling of the original 
bands. Nanotubes close to half filling are deemed to be Luttinger liquids down to milli-Kelvin temperatures 
[5] [6]. In this letter we use the Hubbard Hamiltonian H on the honeycomb lattice to represent the valence 
bands of carbon single- wall (N, N) nanotubes and propose a symmetry-driven configuration interaction pairing 
mechanism which works away from half-filling. We present analytic expressions for the effective interaction 
and obtain the binding energy for (N, N) armchair nanotubes; in the case N = 1 we verify these analytic 
results by exact diagonalization and get high- precision agreement. Starting from the undoped system we find 
that the pair binding energy grows as the number of electrons per C atom increases. Furthermore, we obtain 
stronger binding in nanotubes than in graphite sheets and this suggests a higher critical temperature for the 
former. This is also supported by the measurements of a T c ~ 15 K in the 4 Angstrom single-wall nanotube 
(SWNT) by Tang et al. [7]. 

Using standard notation, the full Hamiltonian reads 

H = H +W = t'J2'J2 ( c Lc r /, CT + h.c) +U^2 nr,]fi rA , (1) 

(r,r') cr r 

where r denotes the honeycomb site, the sum runs over the pairs (r, r') of nearest neighbour carbon atoms and 
t is the hopping parameter. The one-body eigenvalues e (k), (- for the bonding and + for the antibonding 
bands) are readily obtained, and the Fermi line has Ci v symmetry for the nanotubes {Cq v for the graphite 
sheet). Here, we assume that the Fermi level ef lies in the + band. The Hamiltonian in Eq.(l) admits two- 
body singlet eigenstates with no double occupancy and we shall refer to them as W = pairs. The particles 
forming a W — pair have no direct interaction and are the main candidates to achieve bound states in purely 
repulsive Hubbard models already used for the Cuprates [8] [9] [10] . We note that such states are also building 
bricks of the ground state of Hubbard and related models at half filling [11] [12] [13]. 

Recently we obtained [10] [14] a general criterion to get all the W = pairs. We can do that in terms of 
the Optimal Group Q of the Hamiltonian, that we define as a symmetry Group that justifies the degeneracy 
of the single particle energy levels. We may say that an irreducible representation (irrep) r\ is represented in 
the one-body spectrum of H if at least one of the one-body levels belongs to n. Let £ be the set of the irreps 
of Q which are represented in the one-body spectrum of H (£ includes all the irreps when Q is Abelian). Let 
be a two-body eigenstate of the kinetic energy Ho with vanishing z component of the spin. Then, it holds 
the W=0 Theorem: 

77 i £ oWP iv) \ij) =0 (2) 

where is the projection operator on the irrep rj. In other terms, any nonvanishing projection of \ip) on an 
irrep not contained in £ , is an eigenstate of Ho with no double occupancy. The singlet component of this state 
is a W — pair. Conversely, any pair belonging to an irrep represented in the one-body spectrum must have 
positive W expectation value. If a subgroup of Q is available, the => implication still holds. Using the Space 
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Group, we find that in the vanishing quasi-momentum sector the only W = pairs belong to the pseudoscalar 
irrep A^. Let (a, b) denote the basis of the Bravais lattice and u (k, £) the periodic part of the Bloch function 
of quasi-momentum k, with £ = a,b. The pair wavefunction reads [15] 
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with Ri = the origin of the cell where the particle i lies. We can verify by direct inspection that 

ip^ 2 ) 2 (k, Ri, R2) vanishes for Xi = X2, that is the two-body singlet wavefunction vanishes if the particles lie 

on the same annulus of the (N,N) tube. As a consequence Wwl, s , Ri, R2) is an eigenstate of the kinetic 
energy Ho [with eigenvalue 2e(k)] and of the on-site Hubbard repulsion W with vanishing eigenvalue of the 



latter, that is ip\^ 2 } (k, Ri,R2) is a W = pair. Remarkably, ip 
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when the transverse component 



The effective interaction W e g between the particles of a W — pair can be obtained analytically by 



a canonical transformation in the spirit of Ref. [9]. Letting riu (p) denote the non-interacting occupation 



number in band v with wavevector p, we find 
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where X^iO) is the character in 77 of the operation O of C-i v , E is the interacting pair energy and 
[^(ki, k2, k3, k4) is the interaction vertex, with incoming legs k3 and k4 in band + and outgoing ki in 
band + and k 2 legs in band v. The effective Schrodinger equation for the pair reads 



[2e(k) + W F + F{k,E)]a k + ^ W e g (k, k', E) a* = Ea k 

k'£X>/4 



(•») 



where Wf is the first-order self-energy shift and 



[1 - n^> (k + p - q)] [1 - nf (q)] rtf> (p) 
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F(k,£) = -2<5(k-k') 
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is the forward scattering term which does not contains any direct interaction between the particles of the pair. 
Eq.(5) requires a self-consistent calculation of E (since W e s and F are independent). The indices k and k' 
run over 1/4 of the empty part of the FBZ and we denoted such a set of wavevectors as D/4. We show below 
that E = 2ef + Wf + F m i n (kp) + A, with a positive binding energy —A of the W — pair; here -FminCkp) is 
the minimum value of F(k,E) among the kF-wavevectors on the Fermi line. 
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FIG. 1. Comparison between A and A(4) in units of t versus U/t. 



We got a direct verification that pairing actually occours by exact diagonalization for the (1, 1) nanotube 
of length L = 2 (in units of the lattice spacing) and periodic boundary conditions. We define, following Refs. 
[8] [16], 

A(AA + 2)=£(AA + 2)+£(AA)-2£(AA+l), (7) 
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where E(M) is the ground state energy with M electrons (referenced to the electron vacuum); | A (A/" + 2)| is 
one definition of the pairing energy. In previous studies of W = pairing in finite systems we found [8] [10] 
[17] that at least at weak coupling A agrees well with A as obtained by the canonical transformation. For the 
(1, 1) nanotube with N = 2, we can see in Fig.(l) that the agreement between A(4) and A is again very good 
up to U/t m 1. However, we emphasize that A(4) decreases up to a characteristic value of U/t ~ 4 -f- 5, where 
a minimum is reached [15]; at the minimum A(4) ~ —0.018 t. The ratio of the second derivatives with respect 
to U/t at U = was estimated by using best fits and turned out to be 1.00003, while the first derivative 
vanishes. The binding energy for U < t is in the 10~ 3 i range and by analysing the four-body ground-state 
wavefunction we assessed that the symmetry is E\ as predicted by the canonical transformation. Here, E\_ is 
a twice degenerate irrep of the Optimal Group [15] Q that breaks into Ai © Bi in Civ This result encourages 
us to proceed with larger and more physical systems. 
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Table I. Pair binding energy —A (in units of I0~ 3 t) and Average Effective Interaction V (in units of t) for (TV, N) 
nanotubes of various lengths L, as a function of the Fermi energy ep (in units of t). Numerical values were computed 
U/t = 2.5 for illustration. 

We considered supercells of 2 N x L = Nc cells, where L is the length of the (N, N) nanotube in units 
of the lattice spacing. We solved the Cooper-like equation in a virtually exact way for N up to 6 and L 
up to 25, using U/t = 2.5 (which is of the correct order of magnitude for graphite [18] [19]). The canonical 
transformation overestimates A in this range of U/t, but remains qualitatively correct. 

The calculations are performed with the Fermi energy ef varying between 0.8 t and 1.1 t (half filling 
corresponds to ef = 0). As in the (1, 1) cluster, the W = singlets show pairing (see Table I) albeit in general 
in 1 Ai, as expected. The binding energy —A of the pairs decreases monotonically both with the radius and 
the length of the tube. 




<") (b) 
FIG. 2. (a) Results of the canonical transformation approach with U/t = 2.5. — A aS ympt as a function of the Fermi 
energy ef for N = 4 (black boxes), N = 6 (empty triangles) and N = 10 (grey diamonds). The Fermi energy varies 
in the range 0.8 -r 1.1 t. (b) — A asymp t as a function of N for N in the range 6-^36 with ep = t and Average Effective 
Interaction V = 1.5 t. In both figures —Aasympt is in units of t. 

With supercell sizes Nc > 300 numerical calculations become hard. Since we are concerned with the 
asymptotic behaviour for fixed TV and L — + oo and A(N, L) depends on N and L in a complicated way, we 
need a method to make reliable extrapolations of the numerical results. To this end, like in previous work [9], 
[20] we define the Average Effective Interaction V . This is such that setting in Eq.(5) W e s = — jf^, with a 
constant V > for all k and k' in T>/A, one obtains the correct value of A. In other terms, once the binding 
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energy —A(N, L) is known, the constant V must be chosen in such a way that 

V = Nc [2e(k) + F(k)] - [2e F + F min (k F )] - A(N, L) ' (8) 

In Table I we have reported V values; these remain fairly stable around « 1.5-^2 t for A?" > 2 with increasing 
L. Therefore V is largely independent on the Fermi energy and on the radius and this allows us to extrapolate 
to Aasympt^iV) = limz,^oo A(N, L). For N = 4 and N = 6 we use for the average effective interaction V 
the arithmetical mean of the V values reported in Table I for L — 10, 15 and 25; the results are shown in 
Fig. (2. a) together with A asymp t(10) computed with V = 1.5 t. We found that A asympt is strongly dependent 
on the filling at fixed N; the sharp maximum at the optimal doping ef ~ t (which corresponds to a number 
of electrons per graphite atom of 1.25) can be understood in terms of a corresponding peak in the density of 
states. In the optimally doped case — A asymp t(A r ) decreases monotonically as the radius of the tube increases, 
sec Fig. (2.6). The decreasing of the binding energy with N is suggested by recent measurements on nanotubes 
with diameter of few Angstrom [7]. However, in the limit of large N, A asym pt(A'') remains stable around 0.0028 
t and may be interpreted as the binding energy of the W = pair in an optimally doped graphite sheet. 

The paired state we have obtained here is essentially two-dimensional, that is the transverse direction is 
crucial to have a non-Abelian symmetry group and hence W = pairs; the pairing mechanism uses degenerate 
electronic states that exist in 2d and works away from half filling. This opens up the interesting possibility 
that in nanotubes two distinct superconducting order parameters appear in the phase diagram, if it turns out 
that close to half-filling there is another one due to a breakdown of the Luttinger liquid [21]. 

Currently, intercalated graphite and carbon nanotubes superconduct at much lower temperatures than high- 
T c Cuprates and the two kinds of materials are apparently quite different. However, symmetry arguments 
based on the W = theorem tell us that, despite the obvious differences, part of the story must be the same, 
i.e. by a suitable choice of Dirac's characters the on-site Coulomb interaction is utterly turned off. This 
produces the singlet pairing and constrains the ground state spin-orbital symmetry of the interacting system. 
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